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Abstract 

A  wireless  packet  network  is  considered  in  which  each  user  transmits  a  stream  of  packets  to  its  destination. 
The  transmit  power  of  each  user  interferes  with  the  transmission  of  all  other  users.  A  convex  cost  function  of  the 
completion  times  of  the  user  packets  are  minimized  by  optimally  allocating  the  users’  transmission  power  subject 
to  their  respective  power  constraints.  It  is  shown  that,  at  all  ranges  of  SINR,  completion  time  minimization  can 
be  formulated  as  a  convex  optimization  problem  and  hence  can  be  efficiently  solved.  When  channel  knowledge  is 
imperfect,  robust  power  control  is  considered  based  on  the  channel  fading  distribution  subject  to  outage  probability 
constraints.  The  problem  is  shown  to  be  convex  when  the  fading  distribution  is  log-concave  in  exponentiated  channel 
power  gains;  e.g.,  when  each  user  is  under  independent  Rayleigh,  Nakagami,  or  log-normal  fading. 


I.  Introduction 

Power  control  plays  a  crucial  role  in  the  operation  of  a  wireless  network,  as  it  strives  to  provide 
maximum  benefits  to  the  users  within  the  confines  of  available  resources.  In  particular,  in  a  wireless 
network,  each  user’s  transmit  power  interferes  with  the  transmission  of  all  other  users;  therefore,  power 
allocation  has  a  significant  impact  on  the  quality-of-service  (QoS)  experienced  by  the  network  users. 
The  user  benefits  derived  from  a  given  power  allocation  assignment  can  be  characterized  by  different 
performance  metrics.  Traditionally,  a  common  performance  metric  is  a  utility  function  of  the  rates  of  the 
users. 

However,  maximizing  a  concave  utility  function  over  the  feasible  rate  region  in  a  wireless  network  is  not 
necessarily  a  convex  optimization  problem  [l]-[5].  In  the  high  signal  to  interference-plus-noise  (SINR) 
regime,  rate  utility  maximization  can  be  approximately  formulated  as  a  convex  optimization  problem  [1], 
and  hence  it  can  be  efficiently  solved.  For  voice  or  video  applications  that  need  to  maintain  at  least  a 
moderate  minimum  rate,  the  high-SINR  regime  is  often  an  appropriate  assumption.  On  the  other  hand,  in 
wireless  sensor  networks  or  delay-insensitive  data  applications,  a  user  may  wish  to  transmit  at  arbitrarily 
low  rates  (i.e.,  at  low  SINRs)  depending  on  the  channel  conditions.  At  medium  to  low  SINR,  rate  utility 
maximization  are  handled  by  iterative  approximation  methods  [1],  [2].  In  this  paper,  we  study  a  different 
network  performance  metric  that  is  motivated  by  data  applications.  In  particular,  we  consider  the  scenario 
where  each  user  transmits  a  stream  of  packets  to  its  destination,  and  we  wish  to  minimize  a  convex 
cost  function  of  the  user  packet  completion  times.  We  show  that  completion  time  minimization  and  the 
corresponding  optimal  power  allocation  can  be  formulated  as  a  convex  optimization  problem  at  all  ranges 
of  SINR. 

Moreover,  we  consider  imperfect  channel  knowledge  due  to  channel  fading,  and  formulate  the  mini¬ 
mization  as  a  stochastic  programming  problem.  The  channel  gains  are  modeled  as  random  variables:  an 
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outage  event  occurs  when  the  SINR  realization  falls  below  the  transmitter’s  SINR  target.  Robust  power 
control  is  considered  where  each  user  is  subject  to  an  outage  probability  constraint.  We  show  that  for 
a  wide  class  of  commonly  used  channel  fading  distributions,  e.g.,  Rayleigh,  Nakagami,  and  log-normal, 
robust  power  control  can  be  posed  a  convex  optimization  problem. 

Optimal  power  control  in  wireless  networks  is  studied  in  [1];  it  shows  that  maximizing  concave  rate 
utility  functions  can  be  formulated  as  geometric  programming  (GP)  problems  in  the  high-SINR  regime, 
which  are  convex  and  hence  their  solutions  can  be  computed  efficiently.  In  the  medium-  to  low-SINR 
regime,  [6]  describes  an  iterative  approximation  method  to  maximize  a  concave  rate  utility  function 
by  solving  a  series  of  GPs,  and  [2]  proposes  an  algorithm  to  maximize  a  linear  function  of  the  rates  by 
approximating  the  feasible  SINR  region  by  a  series  of  polyblocks.  In  [3],  sufficient  conditions  are  presented 
for  the  convexity  of  the  feasible  QoS  region,  with  optional  constraints  on  the  allocation  of  user  power. 
Utility  maximization  through  joint  optimization  of  adaptive  modulation,  rate  allocation,  and  power  control 
is  investigated  in  [5].  In  [7]  power  control  is  studied  in  frequency-selective  Gaussian  interference  channels. 
Outage  probabilities  corresponding  to  different  fading  distributions  for  network  users  and  interferes  are 
derived  in  [8].  For  interference-limited  wireless  networks,  optimal  power  control  is  considered  in  [9]  under 
Rayleigh  fading  subject  to  outage  probability  constraints. 

The  rest  of  the  paper  is  organized  as  follows.  Section  HD  describes  the  channel  model  and  the  optimization 
framework.  In  Section  [in]  minimizing  a  cost  function  of  the  completion  times  is  posed  as  a  convex 
optimization  problem.  Section  [IV]  explains  how  the  completion  time  region  and  its  corresponding  rate 
region  are  related.  Section  [V]  considers  robust  power  control  against  imperfect  channel  knowledge  subject 
to  outage  probability  constraints.  Conclusions  are  presented  in  Section  IVTl 


II.  System  Model 

A.  Wireless  Channels 

Consider  the  scenario  in  which  M  users  are  communicating  in  a  wireless  packet  network.  Each  user  i 
consists  of  a  transmitter  i  and  a  corresponding  receiver  i,  where  i  =  1, . . . ,  M.  Transmitter  i  wishes  to 
send  a  stream  of  equal-length  packets  to  receiver  i,  where  each  packet  has  L,  bits.  We  assume  the  complex 
baseband  additive  white  Gaussian  noise  (AWGN)  channel  model  between  the  transmitters  and  receivers: 

M 

Y  =  v77“Y,  +  Zu  (1) 

3= 1 

where  V,  is  the  observed  signal  at  receiver  i,  X3  is  the  signal  sent  by  transmitter  j,  H%]  is  the  channel 
power  gain  from  transmitter  j  to  receiver  i,  and  Z,  is  independent  Gaussian  noise  with  power  N%.  We 
consider  coherent  detection:  thus  H%3  may  be  assumed  real  as  detection  is  unaffected  by  phase  change 
in  the  received  signal.  In  the  subsequent  sections,  we  consider  different  channel  knowledge  assumptions 
where  s  may  represent  known  constants  or  random  variables.  Suppose  transmitter  i  has  transmit 
power  constraint  Pi.  When  transmitter  i  transmits  at  a  power  level  of  Pt  <  Pr,  in  the  capacity  limit,  the 
transmission  rate  /?,  (bits  per  channel  use)  achieved  by  user  i  is  given  by: 

Ri  =  log(l  +  Si),  (2) 


where  log  is  base  2,  and  St  is  the  signal  to  interference-plus-noise  ratio  (SINR)  at  receiver  i.  When  the 
interference  from  the  other  transmitters  is  treated  as  noise,  the  SINR  at  receiver  i  is: 


Si 


HaPj 

JV«  +  Ej*  HaPi ' 


(3) 


We  assume  each  user  has  a  sufficient  number  of  packets  backlogged  at  the  transmitter  so  that  it  is  always 
transmitting,  and  we  assume  L,;  is  sufficiently  large  to  allow  transmission  close  to  the  channel  capacity 
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Fig.  1.  An  AWGN  wireless  packet  network  with  M  =  2  users. 


limit.  The  completion  time  (number  of  channel  uses)  of  the  transmission  of  each  of  user  V s  packet  is 
given  by: 

Tt  =  Li/  Ri.  (4) 

For  example,  an  AWGN  wireless  packet  network  with  M  —  2  users  is  illustrated  in  Fig.  [0 


B.  Completion  Time  Cost  Function 

We  consider  the  problem  of  minimizing  a  convex  cost  function  of  the  completion  times  T), . . .  ,TM 
by  optimally  choosing  the  users’  transmission  power  subject  to  the  power  constraints:  0  <  Pi  <  Pi, 
i  =  1, . . . ,  M.  Let  J(T)  be  the  completion  time  cost  function,  where  T  =  [T\  . . .  Tm]  ,  and  other  vectors 
are  denoted  similarly  in  this  paper.  We  assume  J(T)  is  jointly-convex  in  Ti, . . .  ,Tm;  the  convexity 
penalizes  overlong  completion  times.  For  example,  the  following  completion  time  cost  functions  are 
convex  [10]: 

^r(T)  =  T[i]  +  ■  •  •  +  T[r]  (5) 

Jp(  T)  =  (mr  +  •  •  •  +  (: TMy)l/p ,  p>  1.  (6) 

In  ©,  T[ji  denotes  the  Ah  largest  component  of  T.  Thus  the  cost  function  Jr(T)  is  the  sum  of  the  r  longest 
completion  times.  As  special  cases,  r  =  1  represents  the  maximum  completion  time:  maxjTi, . . . ,  Tm}, 
and  r  =  M  represents  the  sum  of  the  completion  times:  Xlfli  G-  hi  ©>  the  cost  function  Jp{ T)  is  the 
lp- norm  of  the  user  completion  times. 


III.  Packet  Completion  Time  Minimization 


A.  Perfect  Channel  Estimation 

We  first  consider  the  scenario  where  the  channel  power  gains  H,fs  can  be  accurately  estimated  and 
they  are  known  by  all  users: 


Hi:j  =  Gij,  i  =  1, . . . ,  M,  j  =  1, . . . ,  M,  (7) 

where  Gy’s  are  known  constants.  In  this  case,  the  minimization  of  the  completion  time  cost  function 
J(T)  can  be  mathematically  formulated  as  the  following  optimization  problem: 


minimize  J(T) 
over  Ti,Ri,Pi 
subject  to 

Ti  >  Li/ Ri 


Ri  <  log  ll  + 

0  <  Pi  <  Pi, 


GuPi 


Nr  +  Y.&  GijPi 


(8) 

(9) 

(10) 


(11) 
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where  i  —  1, . . . ,  M,  and  the  problem  data  Gl3  e  K+,  P?:,  N, ,  Lj  e  M++  are  given.  The  constraint  (fTOl)  in 
the  optimization  problem  is  not  convex  [l]-[3],  [5].  However,  we  show  in  the  next  section  (f8l)— dTIT)  can 
be  transformed  into  a  convex  optimization  problem,  and  hence  its  solution  can  be  efficiently  computed. 


B.  Convex  Optimization 

To  formulate  the  the  completion  time  minimization  problem  given  in  (l8T)-(fTTI)  as  a  convex  optimization 
problem,  we  first  rewrite  the  constraints  (HI),  (fTOl)  as: 


Ti  > 
Si  < 


log(l  +  Si)  ’ 

n  p 

i 


i  =  1 

i  —  1, . . . ,  M. 


Next  we  apply  the  change  of  variables: 

Si  =  in  Si,  Pi  =  in  Pi,  i  =  l,...,M, 

where  in  is  the  natural  logarithm.  The  completion  time  minimization  in  (f8l)— (fill)  then  becomes: 


(12) 

(13) 


(14) 


minimize  J(T)  (15) 

over  Ti,  St,  Pi 

subject  to 

T  >  Li 

log(l  T  exp(iSj)) 

ln|exp(5'i  —  Pi  —  In  Ga  +  In  Nj)  +  ^  exp(Sj  -  P*  +  P,  -  In  Glt  +  In  Gi3 

Pi  <  In  Pi, 

where  i  —  1, . . . ,  M.  Note  that  the  constraint  in  (fTTT)  follows  from  rewriting  (fl3l)  as: 


|  <  0  (17) 

(18) 


SiPr^Ni  +  Y,  SiP^Pfi^Gi,  <  1,  (19) 

jy* 

and  taking  logarithm  on  both  sides  after  applying  (fl4l).  The  change  of  variables  is  similar  to  the  transfor¬ 
mation  techniques  in  geometric  programming  (GP)  problems  [1],  In  particular,  the  log- sum-exp  function 
in  constraint  (fT71)  is  convex  [10].  The  convexity  of  (fl6l)  can  be  verified  from  its  second-order  conditions, 
as  the  right-hand  side  of  (fT6l)  is  twice-differentiable  and  its  second  derivative  is  positive: 


dx 2 


[log(l  +  ex)]  1 


ex  In  2  [2e*  —  ln(l  +  ex)\ 
(1  +  ex )2  [ln(l  +  ex )] 3 


>0, 


(20) 

(21) 


which  follows  from  the  inequality  y  >  ln(l  +  y)  for  y  >  0.  Note  that  the  transformation  in  (fl4l)  does 
impose  a  slight  loss  of  generality  as  we  assume  Pj  ^  0.  Nevertheless,  the  formulation  in  (IT5l)— (fl8T)  is 
otherwise  valid  for  all  ranges  of  SINR,  and  its  solution  can  be  efficiently  computed  by  standard  numerical 
methods  in  convex  optimization  [10].  Note  that  we  may  consider  additional  linear  or  convex  constraints 
on  T,  and  sum  power  constraints  on  subsets  of  P:  they  can  be  readily  incorporated  in  the  optimization 
problem  without  violating  its  convexity. 
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C.  Fading  Channels  and  Power  Adaptation 

In  this  section  we  consider  fading  channels,  i.e.,  the  channel  gains  H^’ s  in  ([[])  experience  random 
variations.  In  particular,  we  assume  the  channel  gains  can  be  characterized  by  a  set  of  s  G  {1, . . . ,  S'} 
discrete  fading  states: 


(G*-1)  with  probability  p\ 

G^  with  probability  ps 


s 

=  1,  (22) 

S=1 


where  H  G  R+xM  =  [H^]  is  the  channel  power  gain  matrix,  and  G  G  M+ xM  =  \GtJ\  is  defined 
similarly.  For  example,  the  discrete  states  may  represent  a  finite  set  of  quantized  channel  estimates.  We 
consider  slow  fading  where  the  duration  of  a  fading  state  is  long  compared  to  the  packet  completion 
times.  We  assume  the  channel  state  s  can  be  accurately  estimated  and  it  is  known  by  all  users. 

We  first  consider  the  case  where  each  user  can  adapt  its  transmission  power  level  according  to  the 
fading  state.  Suppose  user  i  transmits  at  power  level  Pj's>  in  fading  state  s,  subject  to  the  average  power 
constraints: 

s 

E [Pi\±'22p8pW  <PU  i  =  (23) 

S=1 


To  minimize  a  cost  function  of  the  expected  completion  times,  the  optimization  problem  can  be  formulated 
as: 


minimize  J(E[Ti], . . . ,  E [TM]) 
over  E[Tj],  T^s\  S-s\  P^s) 

subject  to 


E[r,]  =  5>sr, 


M 


S— 1 


log  ( 1  +  Sf) 


(s)  p(s ) 

a  r  i 


(«)  d(«) 


V  +  E 

jO) 


o  <p; 

jrP,pp<p„ 


S=1 


(24) 

(25) 

(26) 

(27) 

(28) 
(29) 


where  s  —  1, . . . ,  S]  i  —  1, . . . ,  M.  The  optimization  in  (I24l)-(l29l)  can  then  be  transformed  into  a  convex 
optimization  problem  by  similar  techniques  as  described  in  Section  IIII-BI  Note  that  to  minimize  the 
expected  value  of  the  cost  function,  it  can  be  handled  similarly  by  replacing  the  objective  function  in  (l24l) 
by: 


E[J(T)]  = 


8=1 


(30) 


where  convexity  is  preserved  in  the  nonnegative  weighted  sum  of  convex  functions. 
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In  the  case  where  the  user  cannot  adapt  its  transmission  power  level  to  the  fading  state  (i.e.,  the 
transmitter  is  under  a  short-term  power  constraint),  the  optimization  problem  is  similar  to  (I24l)-(l29l).  but 
with  the  average  power  constraint  in  ([29b  replaced  by  separate  power  constraints  for  each  fading  state: 

,  8  =  1,. ..,S,  i  =  l,...,M.  (31) 

Note  that  under  the  short-term  power  constraints  of  (|3TT).  minimizing  the  expected  cost  function  (l30l) 
decomposes  into  S  independent  optimization  problems:  i.e.,  each  of  J{T[S\ . . . ,  ) ,  for  s  =  1 ,S, 

can  be  minimized  separately. 


IV.  Relations  to  Rate  Region 

In  general,  in  a  wireless  network  as  defined  in  ©-©,  minimizing  a  convex  cost  function  J(T)  of  the 
completion  times  Ti, ,  TM  is  not  equivalent  to  maximizing  a  concave  utility  function  U (R)  of  the  rates 
Ri, ... ,  Rm ■  In  particular,  maximizing  f/(R.)  over  the  rate  region  is  in  general  non-convex  [3]:  at  high 
SINR  it  can  be  approximately  formulated  as  a  GR  and  in  the  medium-  to  low-SINR  regime  there  are 
iterative  approximation  methods  [1],  [2].  Suppose  the  cost  function  J+(T)  is  convex  and  nondecreasing  in 
each  argument  7),  then  completion  time  minimization  is  a  special  case  of  rate  utility  maximization  where 
the  optimization  problem  can  be  formulated  as  convex.  To  see  that  minimizing  J+(T)  can  be  posed  as  a 
rate  utility  maximization  problem,  we  define  the  corresponding  rate  utility  function: 

£/r(R)  4  -  J+(Li/Ri,  ■  ■  ■ ,  Lm/Rm).  (32) 


Note  that  minimizing  J+(T)  is  equivalent  to  maximizing  f7r(R),  and  the  utility  function  LV(R)  is  concave 
in  R  as  prescribed  by  the  convexity  composition  rules  [10]. 

Nevertheless,  in  the  converse,  some  rate  utility  maximization  problems  can  be  formulated  as  minimizing 
convex  functions  of  the  completion  times.  Suppose  we  minimize  a  positively  weighted  sum  of  the 
completion  times: 

J„,(T)  =  w \T\  H - b  wmTm,  w  >z  0,  lTw  =  1,  (33) 


where  >z  denotes  component-wise  inequality,  then  it  is  equivalent  to  maximizing: 


Ud{  R) 


w\ 

Rl 


w 


M 


R 


M 


(34) 


where  Ud( R)  is  the  utility  function  that  corresponds  to  minimum  potential  delay  fairness  [11],  with 
w\  =  u',L,,  i  =  1, . . . ,  M.  In  addition,  minimizing  (T)  in  (1331)  is  also  equivalent  to  maximizing  the 
weighted  harmonic  mean  of  the  rates: 


Uh(R) 


w\ 

i+-+ 


w'm\ 

Pm  J 


(35) 


Note  that  by  applying  Jensen’s  inequality  on  the  convex  function  1/a;  for  x  G  M++,  we  have 


1 

w\R,\  H - b  w'mRm 


w\  w'M 

<  —  +  •••  +  —. 
R\  Rm 


(36) 


which  implies: 


w[R\  +  •  •  •  +  w'mRm  2  t4(R)- 


(37) 


Hence  maximizing  (4(R)  provides  a  lower  bound  to  max  u/  R\  +  •  •  •  +  w'mRm,  which  represents  a 
weighted  throughput  of  the  wireless  network.  In  particular,  the  bound  is  tight  when  R\  =  ■  •  ■  =  Rm  as 
equality  is  achieved  in  (1361):  therefore,  maximizing  the  minimum  rate  in  (1381)  below  can  be  formulated  as 
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Fig.  2.  Completion  time  region  (Pi  =  P2  =  OdB). 


a  minimization  of  the  convex  cost  function  Jx( T)  as  given  in  (l39l).  which  corresponds  to  the  maximum 
completion  time: 


f/n(R)  =  min{f?i, . . . ,  RM}  (38) 

Jx( T)  =  max{Ti, . . . ,  TM}.  (39) 


Moreover,  the  entire  rate  region  achievable  under  (fT0l)-(fTTl)  can  be  characterized  in  terms  of  the 
corresponding  completion  time  region.  Specifically,  the  completion  time  region  as  characterized  in  (fl6l)- 
(fl8l)  is  convex,  and  its  boundary  are  given  by  the  minimizer  of  JW(T)  in  (1331)  over  all  w  A  0,  lTw  =  1.  In 
turn,  from  the  monotonicity  of  ©,  each  minimal  completion  time  vector  T*  =  argmin  JW(T)  corresponds 
to  a  maximal  rate  vector  R*  on  the  boundary  on  the  rate  region  (ITOl)— dTTl).  with  R*  =  L%/T* ,  for 
i  —  1, . . . ,  M.  As  a  numerical  example,  we  consider  the  following  2-user  AWGN  wireless  packet  network: 


G 

P 


0.42 

0.63 


0.89 

0.15 


dB, 


L 

N 


10 

10 

O' 

0 


dB, 


(40) 

(41) 


with  maximum  completion  time  constraints:  T,  <  100Lj.  The  completion  time  region  and  its  corresponding 
rate  region  are  shown  in  Fig.  [2]  and  Fig.  [3]  respectively.  Note  that  the  power  constraints  as  given  in  (HTI) 
belong  to  the  low  SINR  regime,  where  the  high-SINR  GP  approximation  does  not  readily  apply.  The 
completion  time  region  is  convex  in  Fig.  |2l  however,  note  that  its  rate  region  counterpart  is  non-convex 
as  can  be  observed  in  Fig.  [3} 


V.  Robust  Power  Control 

A.  Outage  Probability  Constraints 

In  Section  [nil  we  assume  that  the  channel  power  gains  //,? ’s  can  be  accurately  estimated.  However, 
in  a  fading  environment  where  the  channel  estimates  are  updated  not  as  fast  as  the  channels  vary,  the 
transmitters  may  not  know  the  H1]A  perfectly.  In  this  section,  we  consider  the  scenario  where  the  channel 
gains  are  modeled  as  random  variables: 

H  =  W  4  [Wij]  ,  (42) 

where  the  users  know  the  joint  probability  distribution  of  VF^’s  but  do  not  know  their  realization.  As  in 
Section  UlI-Cl  we  assume  the  duration  of  a  fading  state  is  long  compared  to  the  packet  completion  times. 


0.22 


Fig.  3.  Rate  region  (Pi  =  P2  =  OdB). 


Based  on  the  channel  distribution,  each  user  i  chooses  a  target  SINR  S).  Should  the  realized  channel 
SINR  fall  below  the  target  Si,  the  receiver  cannot  decode  the  transmitter’s  message,  and  it  results  in 
an  outage  event.  To  ensure  the  network  operates  with  an  acceptable  level  of  reliability,  we  extend  the 
completion  time  minimization  framework  in  Section  HIT]  to  additionally  consider  constraints  imposed  on 
the  permissible  probability  of  outage.  Specifically,  we  minimize  the  completion  time  cost  function  subject 
to  a  set  of  outage  probability  constraints:  q,  e  [0, 1],  i  —  1, . . . ,  M,  where  we  stipulate  that  the  probability 
of  user  i’s  transmission  in  outage  not  exceed  g*. 

Incorporating  the  outage  probability  constraints  qf  s,  the  minimization  of  the  completion  time  cost 
function  is  described  by  the  following  stochastic  programming  [12]  problem: 


minimize  J(T) 

(43) 

over  Th  Sh  Pi 

subject  to 

Li 

p.  "»  * 

(44) 

log(l  +  St) 

Pr(  ^Wii  <  S,\  <  qt 

(45) 

0  <  Pi  <  Pi, 

(46) 

where  i  —  1, ,  M.  We  show  that  the  minimization  in  (I43l)-(l46l)  can  be  posed  as  a  convex  optimization 
problem  for  a  wide  class  of  channel  fading  distributions  commonly  considered  in  wireless  communications. 

B.  Robust  Convex  Optimization 

In  terms  of  the  transformed  variables  in  (fl4l).  we  first  define  the  reliability  function  as: 

P)  =  Prjuser  i  not  in  outage}  (47) 


where  (l48l)  follows  from  rearranging  (1451)  with  the  additional  change  of  variables: 

Wij  =  In  Wij,  i  —  1, . . . ,  M,  j  —  1, . . . ,  M. 


(49) 
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Next,  we  characterize  the  reliability  function  in  terms  of  the  channel  distribution.  Let  W,  be  an  M- 
component  non-negative  random  vector  that  corresponds  to  the  zth  row  of  the  channel  matrix  W : 


[Wa  ...  Wm]T, 

(50) 

and  w i  G  K+7  a  realization  of  W,  .  Under  (l49l).  the  transformed  vectors  W,,  w,  are  defined  similarly.  Let 
/w, (w,)  denote  the  joint  probability  distribution  function  (pdf)  of  W,.  Theorem  |Tj  below  describes  the 
sufficient  condition  that  establishs  the  log-concavity  of  $*(,%,  P),  under  which  (I43l)-(l46l)  can  be  posed  as 
the  following  convex  optimization  problem: 

minimize  J(T) 

(51) 

over  T,  ,  ,5),  P, 

subject  to 

^  u 

‘  -  log(l  +  exp(S'j)) 

(52) 

In  P)  >  ln(l  -  qf) 

(53) 

Pi  <  In  Pi: 

(54) 

where  i  —  1, . . . ,  M.  Let  explw,)  denote  component-wise  exponentiation  of  the  vector  w,. 
Theorem  1:  The  reliability  function  <f>i(S'j,P)  is  log-concave  in  Si,  P  if  ./'w,  (expfw,;))  is 

log-concave 

in  w j. 

The  proof  is  given  in  the  appendix.  For  example,  the  above  condition  is  satisfied  when  each  user  expe¬ 
riences  independent  fading  distributed  as  [13]:  i)  Rayleigh  (richly  scattered  environments),  ii)  Nakagami 
(significant  line-of-sight  propagation),  or  iii)  log-normal  (shadowing  due  to  signal  attenuation).  Therefore, 
in  all  these  cases,  completion  time  minimization  subject  to  outage  probability  constraints  can  be  formulated 
as  convex  optimization  problem  (|5TI)-(|541). 

VI.  Conclusions 

In  this  paper  we  consider  minimizing  a  convex  function  of  the  completion  times  of  user  packets  by 
optimally  allocating  transmission  power  in  a  wireless  network.  We  first  focus  on  the  scenario  where  the 
channel  gains  can  be  estimated  accurately  and  are  known  by  all  users.  We  show  that  completion  time 
minimization  can  be  formulated  as  a  convex  optimization  problem,  and  hence  the  corresponding  optimal 
power  allocation  can  be  efficiently  computed.  The  optimization  formulation  is  valid  for  all  ranges  of  SINR, 
which  is  especially  pertinent  for  wireless  sensor  networks  or  delay-insensitive  data  applications  where  the 
users  may  transmit  at  moderate  or  low  SINRs.  Under  fading  channels  with  transmission  power  adaptation 
across  fading  states,  an  average  power  constraint  can  be  incorporated  into  the  optimization  problem.  We 
show  that  completion  time  minimization  is  a  special  case  of  rate  utility  maximization  for  which  the 
optimization  problem  can  be  posed  as  convex.  In  particular,  in  a  wireless  network,  although  the  feasible 
rate  region  is  non-convex,  the  corresponding  completion  time  region  is  shown  to  be  convex.  Finally,  we 
consider  robust  power  control  under  imperfect  channel  knowledge  in  fading  channels.  Completion  times 
are  minimized  subject  to  outage  probability  constraints  over  the  fading  distribution,  and  we  show  that  for 
a  wide  class  of  commonly  used  fading  distributions,  e.g.,  Rayleigh,  Nakagami,  and  log-normal,  robust 
power  control  can  be  posed  as  a  convex  optimization  problem. 

Appendix 

Proof  of  Theorem  [7] 

We  first  introduce,  in  terms  of  the  transformed  variables  in  (l49l).  the  notation  of  W  ,  G  MM_1 
representing  the  interfering  channel  random  vector: 

w_,  4  [wa  . . .  wii+1  . . .  wiM\ T , 


(55) 
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and  w_i  is  a  realization  of  W_*.  The  pdf  of  W_;  is  given  by  the  marginal: 

/w_4(w-i)  =  /  /wjw,;)  dWlt.  (56) 

JWu 

Next,  conditioning  on  W_j  =  w the  reliability  function  in  (l48l)  is  given  by: 

HSi,  P)  =  f  MSi,  P,  w dw_,.  (57) 

)w_i 

For  notational  convenience,  we  write  P,  w_j)  as  a  composition  of: 

0i(£i,P,w_*)  =  F^..(^(^,P,w_i))  /w_.(w_i)  (58) 

Fwui'tiii)  ~  1  _  Fyy..(wn)  (59) 

&(<%,  p,  w_j)  =  ln|exp(5'i  -  P*  +  In  TV*)  +  ^  exp(S'i  -  P*  +  Pj  +  %)  j,  (60) 


where  F^.fwu)  in  (l59l)  is  the  cumulative  distribution  function  (cdf)  of  Wu,  and  (wu)  is  referred  to 
as  its  complementary  cdf,  which  is  a  nonincreasing  function  in  wtl. 

Proof:  In  the  construction  of  QfSi,  P)  in  (l57l).  with  the  application  of  Lemma  [I]  below,  log-concavity 
is  preserved  [10],  [12]  in  the  integration  in  (1571).  (l56l):  multiplication  in  (l58l):  complementary  cdf  in  (l59l): 
and  composition  of  a  logarithmically  concave,  nonincreasing  function  with  a  convex  function  in  (1601).  ■ 
Lemma  1:  Under  transformation  (1491).  f^r  (w,)  is  log-concave  in  w(. 

Proof:  Consider  the  logarithm  of  the  pdf  of  transformed  random  vector  W,: 

ln/w,(*i)  =  ln|  exP  ( ^  ./w,  (exp(wj))  \  (61) 

^  3= 1  * 

M 

=  + ln  (exp(^))  >  (62) 

j=i 

where  log-concavity  of  /w.  (exp(wj))  in  w,  follows  from  assumption  from  Theorem  [0  ■ 
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